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ABSTRACT 

The  Middleton  Class  A  narrowband  non-Gaussian  noise  model 
[9-12]  is  examined.  It  is  shown  that  this  noise  model  (which 
is  known  to  fit  closely  a  variety  of  non-Gaussian  noises)  can 
itself  be  closely  approximated  by  a  computationally  much 
simpler  noise  model.  It  is  then  shown  by  numerical  examples 
that,  for  the  problem  of  locally  optimum  detection,  the  sim¬ 
plest  form  of  this  approximation  yields  nearly  optimal 
(asymptotic)  performance.  The  performance  of  other  simple 
suboptimal  threshold  detectors  in  Class  A  noise  is  also  exam¬ 
ined.  Finally,  a  useful  relationship  between  the  Class  A  model 
and  the  e-mixture  model  is  developed. 


L  Introduction 


For  many  communication  problems  the  usual  Gaussian  noise  assump¬ 
tion  is  inadequate.  Often  this  is  due  to  the  occurrence  of  low  probability, 
high  amplitude  "spikes".  This  impulsive  component  of  the  interference 
has  been  found  to  be  significant  in  many  problems.  Examples  include 
atmospheric  noise,  where  lightning  discharges  in  the  vicinity  of  the 
receiver  can  cause  such  spikes  [1-3,13-17],  and  underwater  problems 
such  as  sonar  and  submarine  communication,  where  the  ambient  acousti¬ 
cal  noises  may  include  impulses  due  to  noisy  aquatic  animals  such  as 
snapping  shrimp  [4]  or  impulses  due  to  ice  cracking  in  arctic  regions  [5]. 
In  addition  to  these  natural  non-Gaussian  noise  sources,  there  is  a  great 
variety  of  man-made  non-Gaussian  noise  sources  such  as  automobile  igni¬ 
tions,  neon  lights,  and  other  electronic  devices  [8-11]. 

Various  attempts  have  been  made  to  develop  models  of  non-Gaussian 
noises.  These  models  can  be  divided  into  two  groups:  those  which  are 
empirically  motivated  and  those  which  are  physically  motivated.  Empiri¬ 
cal  models  [13,14,18-28]  are  those  developed  to  fit  collected  data,  often 
with  little  regard  for  the  underlying  physical  mechanisms.  Physical 
models  [7-12,15-17],  on  the  other  hand,  attempt  to  model  these  mechan¬ 
isms  directly. 

Among  the  physical  models  of  non-Gaussian  noise  some  of  the  most 
general  are  those  developed  by  Middleton  [7-12].  Middleton  divides  non- 
Gaussian  noise  into  two  classes,  A  and  8.  (There  has  also  been  considera¬ 
tion  of  a  Class  C  which  contains  noises  which  are  sums  of  Class  A  and 
Class  B  components  [9].)  Class  B  noises  are  broadband,  Le.  those  with 


spectra  broader  than  the  passband  of  the  receiver.  Class  A  noises  are 
narrowband,  i.e.  have  spectra  comparable  to  or  narrower  than  the 
receiver  passband. 

In  this  paper  we  examine  the  Middleton  Class  A  noise  model.  We  show 
that,  in  a  wide  variety  of  cases,  the  first-order  noise  probability  density 
function  (PDF)  can  be  closely  approximated  by  PDF’s  having  considerably 
simpler  form.  We  then  consider  the  locally  optimum  (also  called  thres¬ 
hold  or  weak  signal)  detection  problem  in  Class  A  noise.  We  show  that  the 
detector  which  is  locally  optimum  for  the  above  PDF  approximation  per¬ 
forms  very  well  for  the  original  Class  A  model.  We  also  examine  the  per¬ 
formance  of  other  even  simpler  suboptimum  threshold  detectors. 
Finally,  in  Section  V,  a  relationship  is  developed  between  the  Class  A  noise 
model  and  the  e-mixture  model  The  e-mixture  model  is  a  highly  tract¬ 
able  empirical  model,  and  through  this  relationship  the  advantages  of  the 
(physical)  Class  A  model  can  be  carried  over  to  the  mixture  model,  e.g., 
physical  motivation  and  direct  calculation  of  basic  parameters. 
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II.  The  Class  A  Noise  liodel 

Middleton  [7-12]  assumes  the  noise  has  the  form  X(t)  +  N(t)  where 
N(t)  is  a  Gaussian  background  component  and 

*(<)  =  £  (0 

i 

where  Uj  denotes  the  3th  received  waveform  from  an  interference  source 
and  U  is  a  random  parameter.  He  then  assumes  that  the  waveform 
receptions  are  Poisson  distributed  in  time  and  shows  that  the  normalized 
(to  unit  variance)  noise  density  /  (x)  can  be  approximated  canonically  by 


/(*)=  £  Kn /(*;«£)  (2) 

m*0 


where  /  (s  ;o*)  is  the  zero-mean  Gaussian  PDF  with  variance  o*.  The  vari¬ 
ance  of  the  m1*  density  is  given  by 


*  -  rn/A  +  r 
m  -  !  +I> 


and  the  coefficient  is  given  by 


Kn  =« 


(3) 

(4) 


where  A  and  P  are  the  two  basic  parameters  of  the  model.  The  first 
parameter,  A,  is  called  the  "overlap  index"  and  is  defined  by  A  =  vT 
where  v  is  the  rate  of  the  homogeneous  Poisson  process  which  governs 
the  generation  of  the  interfering  waveforms  Uj  and  T  is  the  mean  dura¬ 
tion  of  a  typical  interfering  signal.  The  other  parameter,  T,  is  given  by 
the  ratio  of  the  power  in  the  Gaussian  portion  of  the  interference  to  the 
power  in  the  Poisson  component.  The  corresponding  envelope  distribu- 
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tion  is  given  by 

P(E>E0)  =  J  4  (5) 

m  “0 

for  E(jtO.  Middleton  has  shown  that,  by  adjusting  the  parameters  A  and 
r,  the  density  f  given  in  (2)  can  be  made  to  fit  a  great  variety  of  non- 
Gaussian  noise  densities  [9-12].  Also,  the  parameters  A  and  P  are  physi¬ 
cally  motivated  and  can  be  directly  estimated  (see  [12,9]). 

Unfortunately  the  model  (2)  is  cumbersome.  For  example,  in  [ll], 
Spaulding  and  Middleton  exhibit  the  optimal  nonlinearity  for  detection 
(i.e.  the  likelihood  ratio  /  (x— Sj)//  (x-s0))  and  point  out  that  this  detec¬ 
tor  structure  is  likely  to  be  computationally  burdensome  and  uneconomi¬ 
cal.  Thus  we  would  like  to  develop  detector  nonlinearities  having  simpler 
structure  but  which  retain  the  desirable  properties  of  the  one  given  in 
[11]. 
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m.  Approximating  the  Class  A  llodel 

For  each  Atel  we  define  the  PDF  f  y  to  be  the  (normalized)  M-term 
truncation  of  the  Class  A  noise  PDF  given  in  (2),  i.e. 

*2  Jk/ (*;(£) 

/j|(*)  =  -  (6) 

£  Kn 
m« 0 

m 

Since  the  Km's  are  positive  and  2  XJ*  =  1,  we  have  that  fy  (given  in 

m  »0 

(6))  converges  pointwise  to  f  (given  in  (2))  as  (In  fact,  it  converges 

uniformly.)  Our  goal  in  this  section  is  to  show  that  fy  is  actually  a  very 
good  approximation  to  /  for  small  values  of  M.  Note  that  Ju  ^  a 
weighted  sum  of  a  finite  number  of  Gaussian  densities,  a  model  of  this 
form  was  used  to  model  non-Gaussian  noise  in  [22],  [27]  and  [28]. 

That  the  Class  A  noise  PDF  given  in  (2)  is  a  good  model  for  nar¬ 
rowband  non-Gaussian  noise  has  been  demonstrated  [9-12]  by  showing 
that  (for  appropriate  choice  of  parameters  A  and  F)  the  envelope  distri¬ 
bution  (5)  closely  fits  the  measured  envelope  distribution  of  various 
Class-A-type  noises  (e  g.  interference  from  powerlines  or  machinery).  In 
Figures  1-4  we  have  plotted  the  M-term  (normalized)  truncation  of  the 
envelope  distribution  (5)  for  M=l,2,3,...  We  see  that,  in  each  case,  only 
two  or  three  terms  are  necessary  in  order  for  the  truncated  envelope  dis¬ 
tribution  be  indistinguishable  from  the  true  one.  (We  note  that  a  similar 
observation  has  been  made  by  Berry  [34]  concerning  the  instantaneous 
power  density.)  The  parameters  (A  =0.35,  F=0.0005)  for  Figure  1  are 
used  by  Middleton  [6-12]  to  fit  "interference  (probably)  from  nearby 


powerline,  produced  by  some  kind  of  equipment  fed  by  line"  [0].  Figure 
2’s  parameters  (4=0.0001,  r=50)  fit  data  from  ore-crushing  machinery 
[0-12].  The  parameters  of  Figure  3  (4=0.1,  r=0.00l)  and  of  Figure  4 
(4=0.1,  r=0.1)  are  from  a  range  of  typical  values  [9-12]. 

Figures  1-4  give  rise  to  the  hope  that  a  detector  designed  to  be 
optimal  assuming  the  noise  PDF  is  fy  (with  M  small)  would  perform  well 
when  the  actual  PDF  of  the  noise  is  /.  We  consider  here  the  case  of 
'locally  optimum"  detection  (i.e.  small  signal,  large  time-bandwidth  pro¬ 
duct).  The  details  of  locally  optimum  detection  have  been  presented  in 
many  places  [30-33, 11].  We  will  only  state  the  needed  results. 

Under  mild  regularity  conditions  the  (asymptotic)  performance  (or 
processing  gain)  achievable  using  a  given  detector  nonlinearity  g{x)  when 
the  i.i.d.  noise  has  PDF  h(x)  is  given  by  the  efficacy  functional 

f  iz 


V  (g.h)  = 


f  g(x)h'(x)dx\ 


(7) 


f  g*(x)h(x)dx 


For  a  given  noise  PDF  h(x)  the  locally  optimum  detector  nonlinearity 
g*(x)  is  the  solution  to 

tj  #(fc)  =  max  t )(g  ,h) 


and  is  given  by 


0#(*) 


h(x) 


(8) 


If  we  let  h=/jf  in  (8)  we  obtain  the  locally  optimum  detector  nonlinearity 


-7- 


9m(x)  for  Su  Since  fu  is  sjrmmetric,  gy(x)  is  antisymmetric  (Le. 
9u(~x)s~9u (*))•  In  Figures  5-8,  gZ(x)  is  plotted  for  the  examples  of  Fig¬ 
ures  1-4,  respectively,  for  *>0  and  H=2,3,4,...  (for  11=1,  gi{x)  is  a 
straight  line  with  slope  l/of  =(l+r)/Iy).  We  see  that  in  each  case  g£(x) 
closely  approximates  the  locally  optimum  nonlinearity  for  / . 

Returning  to  the  general  case,  we  have  that  for  any  PDF  h(x)  the  per¬ 
formance  of  its  locally  optimum  nonlinearity  g*(x)  is  given  by 

=T,(pVi)  =  /  h(x)dx.  (9) 

i i*(h)  is  also  known  as  Fisher’s  Information  for  h.  An  interesting  (and 
well  known)  fact  about  the  function  r\  *(h)  is  that  it  is  minimized  by  the 
Gaussian  PDF.  In  fact,  the  locally  optimum  detector  nonlinearity  for 
Gaussian  noise  is  the  linear  detector  (p(x)=x/t72)  which  has  performance 
equal  to  unity  for  all  noise  PDF’s. 

Since  the  Class  A  noise  PDF  /  in  (2)  is  highly  non-Gaussian  we  will 
often  have  rj*{f  )»1.  We  have  seen  that  for  very  small  values  of  M,  fy 
and  gu(x)  closely  approximate  /  and  its  locally  optimum  nonlinearity 
g/(x).  So  there  is  reason  to  believe  that  Tj[gZ,f)  should  be  close  to 
r)  ♦(/ ).  Table  1  bears  this  out.  For  each  of  the  examples  which  we  have 
been  considering  (which  cover  a  fairly  wide  range  of  realistic  values  of  the 
parameters)  we  see  that  the  processing  gain  achievable  using  gi  is 
extremely  close  to  that  achievable  using  g}. 
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IV.  Other  Simple  Suboptimum  Nonlinearities 

In  Section  HI  we  plotted  locally  optimum  (and  suboptimum)  detector 
nonlinearities  for  the  Class  A  model  using  logarithmic  scales.  In  Figure  9 
the  locally  optimum  nonlinearity  for  the  example  used  in  Figures  1  and  5 
(4=0.35,  r=0.0005)  is  replotted  using  linear  scales.  The  vertical  dashed 
lines  at  x  =  0.06  and  x  =0.10  divide  the  x-axis  into  three  regions 
Si  =  ||x|«0.06j,  Sz  =  {0.06<!x  1^0.10},  and  S3  =  {|*  I >0. 10}  which  are 
the  regions  where  gj(x)  is  approximately  linear  (Si),  returning  to  zero 
(S2),  and  approximately  zero  (S3).  Evaluating  the  probability  under  /  of 
each  of  these  regions  (or,  more  intuitively  the  fraction  of  the  data  we 
should  expect  to  fall  in  each  region)  we  have  that 
Pr(Si)  »  0.71,  FVfSz)  »  0.01.  and  Pr[Ss)  w  0.28.  Thus  we  see  that  all  but 
about  1%  of  the  time  the  data  will  fall  in  the  approximately  linear  region 
or  the  approximately  zero  region.  This  leads  us  to  believe  that  gj  can  be 
closely  approximated  by  a  blanker  gi  (also  called  a  hole  puncher)  which 
is  shown  in  Figure  10a.  For  comparison  we  have  also  examined  the  per¬ 
formance  of  a  soft  limiter  gk  (or  clipper)  and  a  hard  limiter  gh  (or  sign 
detector)  which  are  shown  in  Figures  10b  and  10c,  respectively. 

In  Table  2  we  have  given  the  processing  gain  achievable  using  the 
locally  optimum  nonlinearity  gj,  the  blanker  g»,  the  soft  limiter  go.,  and 
the  hard  limiter  gut  (Note  that  the  stars  on  ga  and  get  indicate  that  the 
optimal  value  of  c  is  used).  We  have  included  each  of  the  examples  given 
earlier  as  well  as  two  others.  In  each  case  we  see  that  the  blanker  is 
nearly  optimal  while  the  soft  limiter  and  the  hard  limiter  have  substan¬ 
tially  less  than  optimal  performance.  The  one  exception  to  this  is  the  last 


example  (4=1.0,  F=0. 1)  where  the  soft  limiter  is  nearly  optimal  and  even 
the  hard  limiter  out  performs  the  blanker.  Not  surprisingly  the  locally 
optimum  nonlinearity  for  this  case  (shown  in  Figure  11)  is  more  closely 
approximated  by  a  soft  limiter  than  a  blanker.  We  must  stress  though 
that,  based  on  our  experience,  this  seems  to  be  an  unusual  case.  In  fact 
Table  2  is  quite  representative  of  our  findings  in  general. 

Another  issue  of  importance  when  considering  various  detector  non- 
linearities  is  that  of  robustness  or  sensitivity.  Since  the  blanker  and  soft 
limiter  each  only  depend  on  one  parameter  (the  "cut-off"  parameter  c),  it 
is  fairly  straightforward  to  examine  their  robustness.  (Note  that  a  hard 
limiter  does  not  depend  on  such  a  parameter.)  In  Figure  12  we  plotted 
the  processing  gain  achievable  using  gi  and  gSi  versus  the  cut-off  param¬ 
eter  c  for  the  example  (4=0.35,  r=0.0005)  considered  in  Figures  1,  5,  and 
9.  It  would  seem  from  the  smoothness  and  flatness  of  the  curves  in  Fig¬ 
ure  12  near  their  respective  maxima  that  both  the  soft  limiter  and  the 
blanker  are  quite  insensitive  to  variations  in  the  cut-off  parameter. 

On  the  other  hand,  in  Figure  13,  the  same  two  curves  are  plotted 
using  a  different  scale  on  the  abscissa.  This  new  scale  is  not  the  cut-off 
parameter  c  but  the  probability  of  the  set  (— csixsc  j  under  the  Class  A 
PDF  /(*).  As  mentioned  above,  this  can  be  thought  of  as  the  fraction  of 
the  data  which  we  can  expect  to  fall  in  the  linear  region  of  the  detector 
nonlinearity  (cf.  5j  in  the  first  paragraph  of  this  section).  Since  any  esti¬ 
mate  of  c  ♦,  the  optimal  value  of  the  cut-off  parameter,  would  presumably 
come  from  some  version  of  an  empirical  PDF  (see  [12]),  Figure  13  is  likely 
to  be  a  more  reasonable  way  than  Figure  12  to  examine  the  sensitivity  of 
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the  blanker  and  soft  limiter  to  uncertainties  in  estimating  c*  Results 
similar  to  Figure  13  have  also  been  obtained  for  the  Middleton  Class  B 
(broadband)  noise  model  by  Ingram  and  Houle  [35]. 

It  is  not  unreasonable  at  first  thought  to  assume  that  this  change  in 
scale  would  cause  little  change  in  the  relative  smoothness  and  flatness  of 
the  two  curves.  In  fact  Figure  13  shows  quite  strikingly  that  the  blanker 
is  very  sensitive  while  the  soft  limiter  is  very  insensitive  near  their 
respective  maxima.  This  example  is  again  quite  representative  of  a  wide 
variety  of  other  cases. 
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V.  The  E -Mixture  llodel 

In  this  section  we  consider  the  relationship  between  the  e-mixture 
model  (defined  below)  and  Middleton's  Class  A  model.  In  the  Introduction 
we  mentioned  two  broad  classes  of  models  for  non-Gaussian  noise:  physi¬ 
cal  models  (such  as  the  Middleton  models)  and  empirical  models.  One 
commonly  used  empirical  model  is  the  e-mixture  (or  e-contaminated) 
model  in  which  the  first-order  noise  PDF  has  the  form 

/*(*)  =  (l-c)/o(*)  +  c/i(*)  (10) 

where  ee[0,l]  and  f  q  and  /  j  are  PDF’s.  The  PDF  /o  is  usually  taken  to  be 
a  Gaussian  PDF  representing  background  (i.e.  nonimpulsive)  noise  such 
as  receiver  front-end  noise.  Among  the  choices  for  the  "contaminating" 
PDF  /  x  are  various  "heavy-tailed"  densities  such  as  the  Laplacian  or  dou¬ 
ble  exponential  [‘0,19].  Others  have  chosen  f  x  to  be  also  Gaussian  with 
variance  afl  taken  to  be  many  times  the  variance  of  /0,  ojv  The  ratio 
7*  =  a)/a)0  has  generally  been  taken  to  be  between  1  and  100  [20-26]. 
As  with  other  empirical  models  the  disadvantage  of  the  mixture  model  is 
that  the  parameters  (e  and  7*)  are  not  directly  related  to  the  underlying 
physical  situation  and  hence  are  difficult  to  determine.  The  primary 
advantage  of  the  mixture  model  is  its  analytic  simplicity. 

Ve  saw  in  Section  III  that  for  a  wide  range  of  values  of  the  parame¬ 
ters,  A  and  P,  of  Middleton's  Class  A  noise  model  the  (normalized)  sum  of 
the  first  two  terms  of  the  first-order  noise  PDF  give  a  sufficiently  accu¬ 
rate  approximation  to  the  full  PDF  given  in  (2).  In  these  cases  Su  ls  just 
a  "Gaussian-Gaussian"  mixture  as  in  (10)  with  / o(*)=/ (*;*!). 
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/i(*)=/ (*;*?). 


7* 


(11) 


and 


t 


*1  _  _j4_ 

K0+Kx~  1  +  A 


(12) 


Thus  we  see  that,  when  a  Gaussian-Gaussian  mixture  is  used  to  model 
narrowband  non-Gaussian  noise,  it  may  be  possible  to  obtain  e  and  y2 
from  (11),  (12)  and  the  techniques  already  developed  for  determining  A 
and  r  [12,9]. 

Another  interesting  consequence  of  the  relations  (11)  and  (12)  is  that 
(at  least,  for  Class-A-type  noises)  reasonable  values  for  7*  seem  to  be 
between  100  and  104.  For  example,  if  A=0.l  and  r=0.001  as  in  Figures  3 
and  7  then  7Z=10,001,  and  if  A=0.35  and  r= 0.0005  as  in  Figures  1  and  5 
then  72=5715.  These  values  are  as  much  as  two  orders  of  magnitude 
greater  than  those  previously  used  for  7*  [20-26].  On  the  other  hand  the 
commonly  used  values  of  e,  say  between  0.01  and  0.25,  correspond 
approximately  to  values  of  A  between  0.01  and  1/3. 
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VI.  Summary  and  Conclusions 

In  this  paper  we  have  shown  that  the  Middleton  Class  A  noise  model 
can  often  be  approximated  closely  by  the  (normalized)  sum  of  just  the 
first  few  terms.  In  fact,  in  many  cases,  two  terms  are  sufficient.  This  was 
especially  clear  when  we  looked  at  the  efficacy  of  the  detector  nonlinear* 
ity  which  is  locally  optimum  for  two  terms  of  the  Class  A  model  and  found 
it  comparable  to  the  efficacy  of  the  full  locally  optimum  detector  (see 
Table  1). 

We  have  also  examined  the  performance  of  some  simple  suboptimum 
detector  nonlinearities.  For  most  Class  A  noises  the  blanker  has  nearly 
optimal  performance  while  the  soft  limiter  and  hard  limiter  have 
significantly  lower  performance.  On  the  other  hand  the  performance  of 
the  blanker  seems  to  be  far  more  sensitive  to  errors  in  estimating  the 
optimal  cut-off  parameter. 

Finally,  for  those  cases  where  two  terms  of  the  Class  A  model  are 
enough  we  have  developed  a  relationship  between  the  Class  A  model  and 
the  e-mixture  model.  This  yields  a  way  of  estimating  the  parameters  of 
the  mixture  model  directly  from  the  physically  motivated  parameters  of 
the  Class  A  model. 

It  should  also  be  noted  that  one  of  the  clear  advantages  of  the  Mid¬ 
dleton  Class  A  noise  model  is  that,  since  its  parameters  can  be  estimated 
directly  from  the  data,  any  detector  based  on  this  model  could  easily  be 
implemented  adaptively.  Such  a  situation  would  require  real-time  com¬ 
putation  of  the  detector  nonlinearity.  Here  the  above  approximation 
would  result  in  substantial  computational  savings.  Furthermore,  in  such 
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a  situation  it  would  be  possible  to  determine  M  (the  truncation  parame¬ 
ter)  in  an  adaptive  fashion  as  well.  That  is,  first  determine  A  and  P  and 
then  choose  the  number  of  terms  necessary  to  achieve  the  desired 
degree  of  approximation.  In  fact,  this  approach  (of  determining  M  adap¬ 
tively)  could  actually  be  used  to  obtain  optimal  (not  just  nearly  optimal) 
performance  in  some  situations.  Optimal  performance  could  be  achieved 
using  the  philosophy  of  the  Schwartz  detector  (see  [29]  for  details). 

Finally  we  note  that,  while  we  have  carried  out  our  analysis  only  for 
the  problem  of  locally  optimum  detection,  the  closeness  of  the  envelope 
distribution  approximations  (Figures  5-8)  encourages  us  to  believe  that 
similar  computational  savings  might  be  realized  by  applying  the  approach 
of  this  paper  to  other  non-Gaussian  signal  processing  problems.  Further¬ 
more,  we  re-emphasize  that,  while  our  results  have  primarily  been 
demonstrated  by  some  examples,  these  examples  cover  a  wide  range  of 
realistic  values  for  the  parameters  of  the  model. 
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Figure  8.  Locally  optimum  detector  nonlinearity  for  noise 
PDF/*.  M=2,3,4.  ..  (A=.3.r=.l). 
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figure  9.  Locally  optimum  detector  nonlinearity  for  Class  A  noise  PDF  /  (linear  scales) 
(A =0.35,  r=0.0005), 


Table  1.  fjigi-  /)•  Processing  gain  (efficacy)  achievable  using  nonlinearity  gi 
in  Class  A  noise.  17(04.  /)  =  £*(/)  in  all  four  examples. 


Table  2.  Processing  gain  (efficacy)  achievable  using  optimal  (p*).  blanker  (p£). 
soft  limiter  (gsi).  and  hard  limiter  {gla)  nonlinearities. 
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